Элементарные дифференциальные уравнения и задачи с граничными условиями: от уравнений высшего порядка к системам первого порядка

Переход от дифференциальных уравнений высшего порядка к системам первого порядка означает кардинальную смену перспективы. Вместо отслеживания ускорения одной переменной мы развиваем вектор пространства состояний, вектор пространства состояний который одновременно представляет положение, скорость и более высокие производные. Любое линейное уравнение порядка $n$ можно разложить на связанную систему из $n$ уравнений первого порядка, что позволяет использовать всю мощь матричной алгебры.

1. Метод понижения порядка

Чтобы преобразовать скалярное уравнение порядка $n$: $y^{(n)} = F(t, y, y', \dots, y^{(n-1)})$, мы определяем набор вспомогательных переменных:

$$x_1 = y, x_2 = y', \dots, x_n = y^{(n-1)}$$

Эта замена приводит к векторному уравнению $\mathbf{x}' = \mathbf{f}(t, \mathbf{x})$. Для классического механического осциллятора, описанного уравнением $$mu'' + \gamma u' + ku = F(t)$$, преобразование даёт:

$x_1' = x_2$
$x_2' = -\frac{k}{m}x_1 - \frac{\gamma}{m}x_2 + \frac{1}{m}F(t)$

Пример 1: Преобразование системы пружина-масса

Задача

Движение определённой системы пружина-масса описывается дифференциальным уравнением второго порядка $u'' + \frac{1}{8}u' + u = 0$. Перепишите это уравнение как систему уравнений первого порядка.

Подстановка

Пусть $x_1 = u$ (положение) и $x_2 = u'$ (скорость). Следовательно, $x_1' = x_2$.

Матричная форма

Подставляя в ОДУ: $x_2' + \frac{1}{8}x_2 + x_1 = 0 \Rightarrow x_2' = -x_1 - \frac{1}{8}x_2$.

$$\mathbf{x}' = \begin{pmatrix} 0 & 1 \\ -1 & -1/8 \end{pmatrix} \mathbf{x}$$

2. Связанные физические системы

Хотя понижение порядка — это математическая удобная процедура для одного уравнения, системы уравнений возникают естественно в сложных средах:

Механические системы: Много-массовые системы (например, рисунок 7.1.1) включают связанные силы, где движение одной массы влияет на другую через закон Гука.
Связанные резервуары: Течение жидкости между резервуарами (рисунок 7.1.6) зависит от закона сохранения массы, при котором скорость изменения количества соли в резервуаре 1 зависит от концентрации в резервуаре 2.
Электрические цепи: Используя фундаментальные соотношения $$V = RI, C \frac{dV}{dt} = I, L \frac{dI}{dt} = V$$, мы строим системы, описывающие одновременное изменение напряжения и тока в индуктивностях (L), конденсаторах (C) и резисторах (R).

🎯 Основной принцип

Обрабатывая производные как независимые переменные в векторе, мы преобразуем сложность «изменения скорости изменения» в геометрическое вращение и масштабирование в пространстве состояний.

ВОПРОС 1

Преобразуйте уравнение $u'' + 0.5u' + 2u = 0$ в систему первого порядка. Какая система уравнений является правильной?

$x_1' = x_2, x_2' = -2x_1 - 0.5x_2$

$x_1' = x_2, x_2' = 2x_1 + 0.5x_2$

$x_1' = -0.5x_1, x_2' = -2x_2$

$x_1' = x_2, x_2' = -0.5x_1 - 2x_2$

ВОПРОС 2

Для начальной задачи $u'' + 0.25u' + 4u = 2\cos(3t)$ с $u(0) = 1$, $u'(0) = -2$, какие преобразованные начальные условия для вектора $\mathbf{x}(0)$?

$x_1(0) = 1, x_2(0) = -2$

$x_1(0) = -2, x_2(0) = 1$

$x_1(0) = 1, x_2(0) = 0$

$x_1(0) = 4, x_2(0) = 0.25$

ВОПРОС 3

В системе взаимосвязанных резервуаров (рисунок 7.1.6), если резервуар 1 сливается в резервуар 2, почему система уравнений становится «связанной»?

Потому что скорость изменения в резервуаре 2 зависит от концентрации соли в резервуаре 1.

Потому что постоянная гравитации одинакова для обоих резервуаров.

Потому что общий объём воды должен оставаться точно постоянным.

Потому что только один резервуар имеет внешний источник соли.

ВОПРОС 4

Какой компонент в цепи RLC обычно связан с уравнением первого порядка $L \frac{dI}{dt} = V$?

Катушка индуктивности

Конденсатор

Резистор

Источник напряжения

ВОПРОС 5

Если матрица $\mathbf{A}$, представляющая систему, имеет фундаментальную матрицу $\Psi(t)$, как вычисляется специальная фундаментальная матрица $\Phi(t)$, для которой $\Phi(0) = \mathbf{I}$?

$\Phi(t) = \Psi(t)\Psi(0)^{-1}$

$\Phi(t) = \Psi(t) + \mathbf{I}$

$\Phi(t) = \int \Psi(t) \, dt$

$\Phi(t) = \Psi(0)\Psi(t)^{-1}$

Вызов: анализ собственных значений системы второго порядка

Преобразование и решение систем

Рассмотрим систему $\mathbf{x}' = \begin{pmatrix} 3 & -2 \\ 2 & -2 \end{pmatrix}\mathbf{x}$. Эта матрица может представлять связанную физическую систему, где массы 1 и 2 взаимодействуют.

Вопрос 1

Найдите собственные значения и собственные векторы для матрицы системы.

Решение: Решая $\det(A-\lambda I) = (3-\lambda)(-2-\lambda) + 4 = \lambda^2 - \lambda - 2 = 0$, получаем собственные значения $\lambda_1 = 2$ и $\lambda_2 = -1$.
Для $\lambda_1 = 2$ собственный вектор $\xi_1 = [2, 1]^T$. Для $\lambda_2 = -1$ собственный вектор $\xi_2 = [1, 2]^T$.

Вопрос 2

Определите фундаментальную матрицу $\Phi(t)$, удовлетворяющую условию $\Phi(0) = \mathbf{I}$.

Решение: Сначала составьте общую фундаментальную матрицу $\Psi(t) = \begin{pmatrix} 2e^{2t} & e^{-t} \\ e^{2t} & 2e^{-t} \end{pmatrix}$.
Вычислите $\Psi(0) = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}$ и её обратную $\Psi(0)^{-1} = \frac{1}{3}\begin{pmatrix} 2 & -1 \\ -1 & 2 \end{pmatrix}$.
Тогда $\Phi(t) = \Psi(t)\Psi(0)^{-1} = \frac{1}{3}\begin{pmatrix} 4e^{2t}-e^{-t} & -2e^{2t}+2e^{-t} \\ 2e^{2t}-2e^{-t} & -e^{2t}+4e^{-t} \end{pmatrix}$.

Вопрос 3

Объясните физический смысл связывающих элементов (-2 и 2) в матрице системы.

Решение: Элементы вне диагонали представляют связь. В системе масса-пружина эти значения соответствуют жёсткости пружины $k_2$, общей для двух масс. Если бы они были нулевыми, массы двигались бы независимо (несвязанно).